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ABSTRACT 

Of 

We combine the available astrometric and photometric data for the 1993 
microlensing event MACHO-LMC-5 to measure the mass of the lens, M = 0.097± 
0.016 Mq. This is the most precise direct mass measurement of a single star 
\ other than the Sun. In principle, the measurement error could be reduced as 

low as 10% by improving the trig parallax measurement using, for example, 
the Space Inter ferometry Mission. Further improvements might be possible by 
rereducing the original photometric lightcurve using image subtraction or by 
obtaining new, higher-precision baseline photometry of the source. We show 
that the current data strongly limit scenarios in which the lens is a dark (i.e., 
brown-dwarf) companion to the observed M dwarf rather than being the M dwarf 
itself. These results set the stage for a confrontation between mass estimates of 
the M dwarf obtained from spectroscopic and photometric measurements and a 
mass measurement derived directly from the star's gravitational influence. This 
would be the first such confrontation for any isolated star other than the Sun. 



Subject headings: gravitational lensing — stars: low-mass 
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1. Introduction 

Direct measurements of stellar masses provide an essential foundation for theoretical 
models of stars. Such measurements must be free of model-dependent assumptions about 
the star's internal physics and so can be obtained only from the star's gravitational effects 
on external objects. Almost all stars with directly measured masses are components of 
binaries. However, binary stars may not always evolve as do single stars, and thus direct mass 
measurements of single stars are of prime importance. The Sun's mass has been measured 
using two completely independent methods: first by applying Newton's generalization of 
Kepler's Third Law to the motion its companions (the planets); second from its deflection 
of light from distant stars (Dyson, Eddington, & Davis 1920; Froeschle, Mignard & Arenou 
1997). For all stars other than the Sun, direct mass measurements have, until recently, only 
been possible with the first method. That is, direct mass measurements have been restricted 
to components of binary systems with orbital periods of decades or less. 

All light-deflection methods rely fundamentally on the equation, 

I — — 4G 8.14 mas , lN 

where 6 E is the angular Einstein radius, M is the mass of the lens, and 7r re i is the lens-source 
relative parallax. Refsdal (1964) was the first to propose that the masses of stars could 
be measured from light deflection. His method was purely astrometric. By measuring the 
relative separation of a nearby lens and a more distant source, one could directly determine 
the relative parallax 7r re i, the angular impact parameter f3, and the maximal deflection AO. 
In the simplest (and typical) case, (3 3> E , 

01 = AO, (/?»0 E ). (2) 

Hence, combining equations (1) and (2) leads to a simple expression for the mass. Refs- 
dal (1964) 's method will be carried out for perhaps a dozen nearby stars using the Space 
Interferometry Mission (SIM, Salim & Gould 2000). 

A second method was proposed by Gould (1992) and first carried out by Alcock et 
al. (1995). Here one uses the accelerated platform of the Earth to measure the microlens 
parallax, 

" E = JM> (3) 

and combines this with an independent measurement of E to obtain the mass. The method 
yields the mass and, simultaneously, the lens-source relative parallax 7r re i: 

M = — , 7T rcl = 7T E E . (4) 

KTl E 
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The major problem here is that the number of events that last long enough for the Earth's 
acceleration to significantly affect the lightcurve is quite small, and for only a small minority 
of these events is it possible to measure 9 E . Indeed, there is only one event with a precisely 
measured mass using this technique (An et al. 2002), and that event is a binary. There are, 
in addition, two single-star events to which this method has been applied, although the mass 
measurements are much cruder. For one, OGLE-2003-BLG-238, the mass is measured to 
only a factor of a few (Jiang et al. 2004). The other is MACHO-LMC-5, which is the subject 
of this paper. 

In addition, Ghosh et al. (2004) have shown that the mass of OGLE-2003-BLG-145/MOA- 
2003-BLG-45 could be determined precisely provided that the lens-source relative proper 
motion is measured astrometrically. One variant of this method, proposed by Delplancke et 
al. (2001), is to measure 7r E using the accelerated platform of the Earth and 6 E from the cen- 
troid shift of the source during the event (Miyamoto & Yoshii 1995; Hog et al. 1995; Walker 
1995). This method would be especially applicable for black-hole microlensing candidates 
(Bennett et al. 2002; Mao 2002; Agol et al. 2002), which are expected to have large enough 
#e to allow such measurement from the ground and which are dark (and so ineligible for the 
relative proper-motion method of Ghosh et al. 2004). 

Finally, this approach can be extended to a greater number of shorter-duration events 
by measuring the microlens parallax tte by comparing the photometric events as seen from 
the Earth and a satellite in solar orbit (Refsdal 1966; Gould 1995) and combining this with 
an astrometric measurement of 6 E using a space interferometer (Boden, Shao & Van Buren 
1998; Paczyhski 1998). Of order 200 such mass measurements will be made by SIM (Gould 
& Salim 1999), which will itself be in solar orbit. 

MACHO-LMC-5 is unique in that it has both a measurement of 7Te derived from the 
microlensing lightcurve (Alcock et al. 2001; Gould 2004), to which we refer in this paper 
as a "photometric" quantity, and a completely independently determined full post-event 
astrometric solution, including both 7r rc i and the lens-source relative proper motion n Tel 
(Alcock et al. 2001; Drake, Cook & Keller 2004). Here, 7r E is the vector microlens parallax, 
whose magnitude is ir E and whose direction is that of the lens-source relative motion. 

Thus, MACHO-LMC-5 is of particular interest for several reasons. First, it permits three 
independent tests on the consistency of the measurements. Second, if the measurements pass 
these consistency tests, they can be combined to obtain a more accurate estimate of the mass. 
Third, the very high magnification of the microlensing event A max ~ 80 permits one to place 
very strong constraints on the lens being a close binary rather than a single star. Fourth, 
the photometric and astrometric measurements can be combined to test the hypothesis that 
the resolved star that appears to be the lens in the astrometric images is not in fact the lens 
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but rather is merely a luminous binary companion to a dark substellar object that generated 
the microlensing event. Finally, the successful completion of all these tests would allow a 
direct confrontation between the measured mass of a single star and its mass as predicted 
from photometric and spectroscopic observations. This would be the first such confrontation 
for any isolated star other than the Sun. 

Here we build on the work of Alcock et al. (2001), Gould (2004), and Drake et al. (2004) 
to carry out the above-described tests, in so far as it is possible today, and outline how these 
tests can be further refined in the future. The error in our mass estimate, 17%, is the smallest 
for any direct mass measurement of a single star other than the Sun, and even approaches 
the precision of measurements of M-dwarf masses from binaries (Delfosse et al. 2000). 

2. Brief History 

Alcock et al. (2001) originally measured n vcl by analyzing epoch 1999 Hubble Space 
Telescope (HST) Wide Field and Planetary Camera (WFPC2) images obtained by Alcock 
et al. (2000). To do so, they assumed that the two resolved objects in those images were 
the lens and the source that had been virtually coincident 6.3 years earlier when the event 
occurred in 1993. They checked for consistency between the direction of /x rel so obtained 
and the direction of 7Te (measured in the heliocentric frame) as determined from a microlens 
parallax analysis of the microlensing-event lightcurve. They found that these position angles 
agreed at the 1 a level. By combining the astrometric and lightcurve data, they determined 
both 7r E and 6*e and so, using equation (4), both the mass and distance of the lens. The best 
estimate of the lens mass was substantially below the hydrogen-burning limit (albeit with 
moderately large error bars) and so appeared somewhat inconsistent with the hypothesis 
that the luminous star seen in the images was actually the lens. The inferred, distance, 
D\ ~ 200 pc, was also quite close relative to the photometric distance inferred from the HST 
photometry. Alcock et al. (2001) suggested that more precise astrometric measurements 
with the new Advanced Camera for Surveys (ACS) on HST could determine the lens-source 
relative parallax 7r re i, and so test the distance estimate and hence effectively test the entire 
procedure for inferring both the mass and distance of the lens. 

In the meantime, Gould (2004), building on the work of Smith et al. (2003), discovered 
that microlens parallax measurements of relatively short events (with timescales £e ^ Y r /27r) 
are subject to a four-fold degeneracy composed of two two-fold ambiguities. While one these 
ambiguities (the so-called "constant-acceleration" degeneracy) has only a very small effect 
on the event parameters, the other (the "jerk-parallax" degeneracy) can affect 7Te by quite a 
lot. Gould (2004) showed that the alternate jerk-parallax solution for MACHO-LMC-5 had 
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a somewhat larger mass and a much larger distance, D\ ~ 450 pc. 

Most recently, Drake et al. (2004) have carried out the ACS measurements advocated 
by Alcock et al. (2001), and these have yielded both a more precise measurement of /x rcl 
and a new measurement of 7r re i- Drake et al. (2004) were able to conduct two consistency 
checks. First, they found that their trigonometric parallax measurement was consistent with 
the lens distance derived by Gould (2004) for his alternate (D\ ~ 450 pc) solution. Second, 
they found that the direction vector ^ rcl //i re i determined from astrometry was consistent 
with the direction vector tt e /7t e found by Gould (2004) for this alternate solution. 

In fact, the proper motion derived from the WFPC2 observations proved very accurate 
and agreed to within 0°.2 with the direction from the new ACS measurements (both in 2000 
celestial coordinates). However, Alcock et al. (2001) reported these results only in ecliptic 
coordinates, and they made a 14° degree error when they translated from celestial to ecliptic. 
This transcription error significantly affected the mass and distance estimates of both Alcock 
et al. (2001) and Gould (2004) when they used the direction of proper motion as a constraint 
in their solutions. The correction of this transformation error, by itself, resolves the most 
puzzling aspects of the solutions obtained by Alcock et al. (2001) and Gould (2004). If the 
correct direction had been incorporated into the fits, both solutions would have moved to 
within 1 a of the hydrogen-burning limit. The Gould (2004) solution would have moved 
to ~ 530 pc and so would have been in better agreement with the photometric parallax, 
although the Alcock et al. (2001) solution would have actually moved to an even shorter 
distance, ~ 160 pc. 



2.1. ACS Astrometry 

In this paper, we will make extensive use of the astrometric measurements of Drake et 
al. (2004), sometimes combining them with other real and/or hypothetical data. To do so, 
we must fit to the original astrometric data. We find that when we fit these data alone, 
the results differ very slightly (much less than 1 a) from the results reported by Drake et al. 
(2004). For consistency, we always use the values from our own fits, although the difference 
has no practical impact on any of the derived results. In fitting the Drake et al. (2004) data, 
we exclude the 2002 F814 point as they also did (A. Drake 2004, private communication). 
This point is a significant (3 a) outlier, which contradicts a F606 measurement taken at 
almost exactly the same time. 

Our fit yields the following parameter estimates: relative parallax 7r rc i = 1.780 ± 
0.185 mas; relative proper motion /z re i = 21.381 ± 0.022 mas yr -1 ; proper motion compo- 
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nents /x re i,East = 17.547 ± 0.029 mas yr 1 ; and // re i,North = —12.217 ± 0.022 mas yr 1 ; and 
position angle = 124°.85 ± 0°.08. 

As noted by Drake et al. (2004), the residuals for the WFPC2 point are quite small 
compared to their reported errors, with A% 2 = 0.11 for 2 degrees of freedom (dof). This 
may imply that the errors were overestimated by Alcock et al. (2001). However, since there 
is a 10% probability of having such low residuals by chance, no definite conclusions can be 
drawn regarding a possible overestimation of the error bars. 

3. Sources and Consistency of Data 

In this paper, we will draw together four sources of data to measure the mass of MACHO- 
LMC-5. We first summarize these sources, then discuss a series of tests that we have carried 
out to determine whether they are consistent with each other. Only after these tests are 
successfully concluded do we combine the data. 

3.1. Data Sources 

The primary data set is the original MACHO SoDoPHOT pipeline photometry of the 
event, which occurred in 1993. These data have already been analyzed by Alcock et al. 
(2001) and Gould (2004). They consist of 352 points in the non-standard MACHO red filter 
(hereafter Rm) an d 265 points in the non-standard MACHO blue filter (hereafter Vm)- We 
slightly deviate from previous authors by recursively removing outliers and renormalizing 
the errors so as to enforce x 2 P er degree of freedom (dof) equal to unity in each bandpass 
separately. We repeat this procedure until all 3.5 a outliers are removed. This removes three 
Rm points and one Vm point, all greater than 3.9 a. The next largest deviation is at 3 a, 
but with more than 600 points, such a deviation is consistent with Gaussian statistics and 
so cannot be considered an outlier. The error renormalization factors are 0.79 in Rm and 
0.81 in V M - 

Gould (2004) had argued against blindly applying this renormalization procedure be- 
cause the mass determination is dominated by the relatively small number of points during 
the event, while x 2 1 dof is dominated by the much larger number of baseline points taken over 
several years. However, we find that x 2 /dof is similar for both these subsets and therefore 
proceed with the renormalization described above. Of the four eliminated points, only one 
is during the event, the R M point at JD = 2490015.14. From Figure 1 of Gould (2004), it 
can be seen that this point is a clear outlier with an abnormally large error bar. 
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For isolated faint stars, SoDoPHOT ("Son of DoPHOT") reports only photon noise 
errors, which of course cannot be overestimates. However, in crowded fields, SoDoPHOT 
follows DoPHOT (Schechter, Mateo & Saha 1993) in "padding" the photon noise as it 
subtracts out surrounding stars. This additional padding is correct in some average sense, 
but may be an overestimate or underestimate in individual cases. Hence, for events of 
particular interest, it is worthwhile to investigate these error bars more closely. 

For Gaussian errors, removal of outliers greater than (T max = 3.5 artificially reduces x 2 
by ~ (2/7r) 1 / 2 o" max exp(— of nax /2) ~ 0.6% and so understates the size of the error bars by 
0.3%. This difference has no practical effect on the results reported here. 

The second data source is the astrometric measurement of the lens-source separation 
made by analyzing epoch 1999 Hubble Space Telescope (HST) Wide Field and Planetary 
Camera (WFPC2) images originally obtained by Alcock et al. (2000). Alcock et al. (2001) 
and Gould (2004) have previously combined the proper motion measurement, /Lt rel , derived 
from these data with the above-mentioned lightcurve data to estimate the mass and distance 
of the lens. 

The third data source is the photometric measurements of the source brightness made 
from the same HST WFPC2 images. These data help constrain the lightcurve fit and thus 
tighten the errors on the mass and distance measurements. Alcock et al. (2001) made use of 
these measurements in their constrained fit, but Gould (2004) did not. 

The final data source is the new astrometric measurements made by Drake et al. (2004) 
using ACS. These have yielded both an improved measurement of /x rel and a new parallax 
measurement, 7r re i- Both the WFPC2 and ACS measurements are listed in Table 1 of Drake 
et al. (2004). 

3.2. Consistency Checks 

3.2.1. Source Color and Magnitude 

We first wish to combine the original SoDoPHOT photometry of the event with the flux 
measurement of the source made by Alcock et al. (2001) after the event was over and the 
source was well separated from the lens. The superb resolution of HST virtually ensures 
that all blended light is removed from the source with the possible exception of a wide- 
binary companion to the source, which we discuss below. This HST source photometry 
can be compared with the source flux that is returned as a parameter by the microlensing 
fit. To do so, one must first translate the HST photometry into the MACHO bands. This 
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can be done by directly comparing the flux levels recorded by HST and SoDoPHOT for an 
ensemble of other stars in the field. While each of these is blended in the MACHO images, 
the blending is equally likely to contaminate the object or the "sky" determination. So it 
should introduce scatter but not a systematic bias. This scatter is smaller for brighter stars, 
but unfortunately the PC field is not big enough to contain many bright stars. Based on 
a comparison of 18 relatively bright stars (and constraining the fits by the two color-color 
slopes, V — V M oc — 0.20(Vm — Rm), R — Rm °c +0.18(Vm — Rm) reported by Alcock et al. 
1999), we find that the HST data imply 

f sy = 23.27 ± 1.34, f BjR = 23.84 ±1.23, ^ = 1.024 ± 0.016 (HST), (5) 

Js,V 

where the errors and covariances are derived by enforcing x 2 /dof in the fit. The error in 
the color calibration is substantially smaller than the errors in the flux calibrations because 
the latter are highly correlated, with correlation coefficient p = 0.962. These results may be 
compared to the source flux levels derived from the fit to the SoDoPHOT data alone, 

f sy = 28.97±3.99, f SiR = 29.50±4.16, ^ = 1.018±0.018 (SoDoPHOT lightcurve). 

Js,V 

(6) 

The two determinations differ by ~ 1.3 o in each of the two (highly correlated) bands sepa- 
rately and by ~ 0.3 a in the color. Since the two photometric measurements are consistent, 
they can be combined. We then find that the x 2 minimum increases by 2.2 for two additional 
dof, confirming the consistency of the two pairs of measurements. 

The one possible caveat is that the microlensing fit gives the flux of the source that was 
magnified during the event while the HST measurement gives all the flux from stars within 
about 100 mas of the source center, corresponding to about 5000 AU at the distance of the 
Large Magellanic Cloud (LMC) (Alves 2004). Significant sources of light that lay beyond 
50 mas would have shown up in subsequent ACS images. Now, as we will show below, the 
angular Einstein radius is #e ~ 1 nias. Any significant light source within about (1/3)#e 
of the source would have betrayed itself during the event. This still leaves the possibility 
that the source has a binary companion between 17 and 2500 AU. While we cannot rule 
out such a possibility, there are several lines of argument against it. First, the discrepancy 
(eqs. [5] and [6]) between the two measurements (though not statistically significant) is of 
the wrong sign to be accounted for by a binary companion. Second, the good agreement in 
the colors shows that any companion must be either of nearly the same color as the source 
or quite faint. If the former, the two stars should also be of roughly the same magnitude, in 
which case one would expect the discrepancy to be much larger than observed (and, again, 
in the opposite direction). Third, to produce ~ 10% or more of the light, the companion 
mass would have to be > 70% of the primary. If LMC binaries are similar to those studied 
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by Duquennoy & Mayor (1991) in the Galactic disk, the fraction of stars with companions 
in the required mass and separation ranges is only ~ 7%. Hence, it is unlikely, though not 
impossible, that such a companion exists and is significantly corrupting the measurement. 

3. 2. 2. Direction of Motion 

We therefore begin by incorporating the WFPC2 photometric measurements into the 
lightcurve fit, taking account of both their errors and covariances. The resulting contour plot 
for the vector parallax 7Te (in the geocentric frame) is shown in Figure 1. This figure should 
be compared to Figure 1 from Gould (2004). Each of the two minima is consistent between 
the two figures at the 1 a level. This is to be expected, since the additional photometric 
data are consistent with those used by Gould (2004). However, the errors are substantially 
smaller, both because the photometric errors have been renormalized by a factor ~ 0.8 and 
because of the additional higher-precision EST baseline photometry of the source. 

We can now ask whether the direction of motion (in the heliocentric frame) implied 
by each of these lightcurve (0i c ) solutions is consistent with the direction of proper motion 
(0 ast ) that we derive from the EST ACS and WFPC2 data of Drake et al. (2004). For the 
southeast solution, which the Drake et al. (2004) parallax measurement demonstrates to be 
the correct one, the comparison yields 

ast = 124°.85 ± 0°.08 lc = 132°.3lfo;|. (7) 

Hence, the photometrically and astrometrically determined directions are consistent at about 
the 1.3 a level. The other (northwest) solution has a direction 4>\c = 136°.5l2°5- However, 
because the x 2 surface deviates strongly from a parabola, the discrepancy with the proper- 
motion data is only at about the 2.8 a level. 

3.2.3. Parallax 

Since the astrometrically determined direction of motion is consistent with the value 
derived from the lightcurve (at least for the southeast solution), we combine the astrometric 
and photometric data. As explained by Alcock et al. (2001) and Gould (2004), this permits a 
full solution for the event, including the mass M of the lens and lens-source relative parallax 
7r re i- We plot the result in the [log M,(M — m) } plane, where (M — m) is the lens distance 
modulus, i.e., corresponding to tt\ = 7i IC \ + 7r s and where we have adopted tt s = 20 /xas for the 
source, which resides in the LMC. Figure 2 shows the resulting likelihood contours. Note 



-10- 



that each set of contours is offset from its respective minimum. The minimum of the short- 
distance (northwest) solution is actually higher by A% 2 = 6, which is a reflection of the mild 
direction discrepancy found for this solution in § 3.2.2. 

Also shown in Figure 2 is the best fit and 1 a error bar for the parallax determina- 
tion based on the astrometric data of Drake et al. (2004), 7r rc i = 1.780 ± 0.185 mas. This 
corresponds to a distance modulus 

(M — m)o — 8.72l°;2i (trig parallax). (8) 

To make an algebraic comparison, we fit the A% 2 < 1 region of the rightward contours to a 
parabola and find, 

log M = -1.023 ± 0.084, (M - m) = 8.683 ± 0.144, p = 0.921, (lc + /x rcl ) (9) 

where p is the correlation coefficient. Hence, the lens distance derived from the lightcurve/proper- 
motion analysis is consistent with the trig parallax at the 1 a level. 

4. Mass, Distance, and Velocity of Lens 

Since the two measurements are consistent, we combine them. The results are shown in 
Figure 3 and can be represented algebraically by, 

logM = -1.013 ±0.073, (M-m)o = 8.702±0.124, p = 0.896. (lc + /x rcl + 7r rcl ) (10) 

These figures correspond to a best-fit mass and distance 

M = 0.097 ± 0.016 M Q , A = 550±30pc. (11) 

This best fit has x 2 — 605.39 compared to x 2 — 601.68 for the lightcurve alone. That is, 
Ax 2 = 3.71 for 5 additional dof. 

Also shown in Figures 2 and 3 are the estimates of the mass and distance of the lens as 
derived by Alcock et al. (2001) from photometric HST data. These are consistent with the 
microlensing measurement at the 1 a level. 

The (U, V, W) velocities of the lens toward the Galactic center, the direction of Galactic 
rotation, and the north Galactic pole are 

U = 43.6 ± 1.9 kms -1 , V = -60.8 ± 8.3 kms -1 , U = 26.6 ± 5.7kms -1 , (12) 

with correlation coefficients puv = —0.90, puw = —0.74, pvw = 0.93, where we have taken 
account of the source motion (van der Marel et al. 2002) and the motion of the Sun relative 
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to the local standard of rest (Dehnen & Binney 1998). The uncertainties are dominated by 
the error in the radial-velocity measurement, v r = 49 ± lOkms^ 1 (Alcock et al. 2001), and 
this fact accounts for the high correlation coefficients. 

Finally, for reference we note that, from equations (1) and (3), these determinations of 
the mass and distance imply 

E = 1.19 ±0.07 mas, vr E = 1.51 ± 0.17, (13) 

where we have taken account of the correlations in determining the errors. 

5. Constraints on Binarity 

An important application of microlensing mass measurements is the opportunity they 
afford to confront theoretical models that attempt to predict the masses of stars from their 
spectroscopic and photometric properties. Crucial to such a comparison is the determination 
that the "star" is in fact a single object and not a close stellar binary or a binary composed of a 
star and a brown dwarf. In the former case, the photometric properties would be a composite 
and in both cases the microlens mass would not be the mass of the star dominating the light. 
It is equally crucial that the luminous star whose visible properties are being measured is 
actually the lens whose mass was measured during the microlensing event, as opposed to a 
luminous companion of a dark object (e.g. a brown dwarf) that generated the microlensing 
event. In this section we investigate how well both of these concerns can be addressed with 
current and/or future data. 

5.1. Limits on Close Binaries 

A close binary approximates a Chang & Refsdal (1979, 1984) lens with sheer 7 = 
[d/(q l l 2 + g -1 / 2 )] 2 , where dO-^ is the angular separation of the two components and q is their 
mass ratio (Dominik 1999; Albrow et al. 2002). These have caustics of full angular width 
47#e, which, if traversed, would give rise to obvious deviations from a point-lens lightcurve. 
For a continuously sampled lightcurve, the best way to avoid the caustic is for the source 
to pass at an angle of 45° relative to the binary axis. In this case, the caustic is just barely 
nicked if Uq = 2 1//2 7, where u is the lens-source separation in units of #e and Uq is its minimum 
value (i.e., impact parameter). For MACHO-LMC-5, the sampling is far from uniform, and 
the tightest simple constraint is obtained from the highest point, which comes almost exactly 
one day after the peak and therefore is at u — Whi g h = 0.029. The largest circle that one can 



- 12 - 



inscribe in the central caustic has radius u ~ 7. From this we derive, 

d 2 , N 

T = (gl/2 + g -l/2)2 < = - 029 ' ( 14 ) 

since otherwise this point would have landed in the caustic and so would have been much 
more magnified than it actually is. Since the lens is much closer than the source, the Einstein 
radius is essentially equal to the projected Einstein radius, = fv = AU/7Te = 0.66 AU. 
Hence the above limit can be expressed in terms of the projected separation between the 
binary components, r±_ = dr E , 

r ± < 0.11(g 1/2 + <T 1/2 ) 2 AU. (15) 

Unless we are very unfortunate to see a widely separated pair projected along the line of 
sight, or unless the companion is of such low mass as to be uninteresting, the putative 
companion would cause the source to move by lOkms -1 or more. This would in principle 
be detectable by spectroscopic measurements. 

While close binaries deviate most sharply from point lenses inside their caustics, they do 
show significant deviations in the surrounding regions as well (see Fig. 1 from Gaudi & Gould 
1997). A detailed accounting of these deviations would strengthen the limit in equation (15) 
but, given the sparse sampling of MACHO-LMC-5, the improvement would most likely be 
modest. 



5.2. Limits on Wide Binaries 



A similar argument places a limit on wide companions, which also give rise to Chang- 
Refsdal caustics. In this case, 7 = qd~ 2 , so 

r± > 3.9g 1/2 AU, (16) 

corresponding to 7.2g 1 / 2 mas. 



5.3. Constraints on the Dark Lens Hypothesis 

As shown in Figures 2 and 3, the best-fit mass lies close to the hydrogen burning limit. 
It is therefore possible in principle that the lens is not actually the red star seen in the 
HST images but rather an invisible brown-dwarf companion to it. To what extent can this 
scenario be constrained by the available data? 
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As pointed out by Drake et al. (2004), even if one relaxes the assumption that the red 
star and the microlensed source were coincident at the time of the event, the remaining 
astrometric measurements "point back" to a relative offset AO very close to zero at the peak 
of the event. Specifically, we find, 

A0n — —1-2 i 7.7 mas A9 E = -2.9 ± 7.5 mas. (17) 

As discussed by Drake et al. (2004), the smallness of these values relative to the errors most 
likely reflects an overestimation of the HST WFPC2 errors by Alcock et al. (2001). However, 
to be conservative, we ignore this possibility. Somewhat stronger constraints can be obtained 
by noting that in this relaxed solution, the parallax error grows from 0.185 mas to 0.291 mas, 
and that this error is fairly strongly correlated with Ad. However, the microlensing analysis 
independently constrains the parallax to be 7r re i = 1.81 ± 0.12 mas, and this constraint can 
be added into the fit. We then find, 

A9 N = -0.8 ± 5.6 mas A9 E = -2.5 ± 6.8 mas. (18) 

These results indicate that, at the 2 a level, the M dwarf must have been within about 13 
mas of the source at the time of the event. On the other hand, from the argument in § 5.2, 
the M dwarf could not have been too close to the source if it were not actually the lens. By 
hypothesis, the putative brown-dwarf lens is below the hydrogen-burning limit while the M 
dwarf is above it, so q > 1. Hence, 

A# = ^>7.2mas. (19) 
A 

Equations (18) and (19) leave only a fairly narrow range of allowed separations. For 
face-on circular orbits and for a total binary mass M tot = 0.2 M , these limits correspond to 
a period range 

18yrs < P < 44yrs (allowed periods). (20) 

Even assuming a mass ratio q = 2, the amplitude of the M dwarf motion would be between 
2.4 mas and 4.3 mas. These amplitudes are quite large relative to the ~ 0.3 mas errors 
achieved for single epoch HST ACS images. Hence, in principle, this scenario could be much 
more tightly constrained by future observations. 

6. Future Confrontations 

From Figure 3, the photometrically derived mass and distance of MACHO-LMC-5 are 
consistent with the values of these properties derived by combining the astrometric and 
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microlensing data. The error bars for both determinations could be improved significantly 
by obtaining additional data and by improving the analysis. 

On the photometric side, both the mass and the luminosity of the M dwarf are inferred 
from its color. Apart from the error in measuring this quantity, these inferences suffer from 
the intrinsic dispersions of mass and luminosity at fixed color. A substantial part of this 
dispersion is due to metallicity. Drake et al. (2004) have argued that the kinematic data 
are consistent with either a disk or a thick-disk star, and therefore a range of metallicities 
of about 1 dex. Hence, spectroscopic determination of the M dwarf's metallicity would go a 
long way toward shrinking the photometry-based mass/distance error bars. 

On the astro metric/microlensing side, there are three paths to improvement. First, of 
course, the distance determination could be improved by a better trigonometric parallax 
measurement. Moreover, because the mass and distance measurements are highly corre- 
lated (see eq. [10]), a more accurate distance would also improve the mass determination. 
Unfortunately, significantly better parallax measurements will not come cheaply. From a 
comparison of equations (8) and (9), the 185 /xas error from the ACS astrometry is about 
50% larger than the distance estimate achieved from microlensing (and the proper-motion 
measurement) alone. A plausible target for a significant improvement would be a 100 /xas, 
or better yet 50 /zas measurement. These would yield mass determinations with fractional 
precisions of 13% and 10%, respectively. 

Note that even if the distance were known exactly, the microlens mass measurement 
error could only be reduced to 7.5%. 

If it were only necessary to consider the statistical errors, such improvements could 
be achieved by multiplying the total length of ACS observations by 3 and 14 respectively. 
However, systematic errors may become important, and the discrepancy between the F606 
and F814 measurements implies that caution is warranted. 

A parallax measurement by SIM might also prove feasible. This seems impossible at 
first sight because the M dwarf has V = 22.7 whereas the magnitude limit of SIM is often 
said to be V < 20. However, what fundamentally limits SIM at faint magnitudes is the 
number of sky photons that enter its 1" radius stop. If we ignore this sky noise for the 
moment, a 50 yuas measurement at V = 22.7 would require an observation of only about 
1 hour. An additional 30- minute observation of the V = 21 source would yield a 30 /xas 
measurement, for a combined error in the relative offset of 60 //as. As in the case of the ACS 
measurements, only two epochs (in two orientations) would be required because it is known 
from the microlensing event that the two stars were virtually coincident in 1993. We find 
that a total of 4 pairs of observations (each pair totaling 90 minutes) would yield a relative 
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parallax error of 42 fxas. The sky is V ~ 22.7 mag arcsec" 2 (so V ~ 21.5 inside the SIM 
stop), which would mean that the observation time for the lens (but not the source) would 
have to be roughly tripled relative to the naive estimates. However, SIM has a very broad 
bandpass and the M dwarf is very red, which may imply that much of the astrometric signal 
will come in well above the sky. Hence, the required duration of the exposures cannot be 
properly estimated until the details of SIMs throughput are worked out in greater detail. 
In any event, it appears that SIM could achieve a substantial improvement in the parallax 
measurement without prohibitive observing time. 

The second potential path would be to obtain better photometry of the source. Recall 
that the HST&ux measurement error was dominated by the problem of aligning the WFPC2 
and SoDoPHOT photometry, which was exacerbated by the small number and faint flux 
levels of stars in the small PC chip. Since the source and lens are now well separated, one 
could image them using the much larger ACS camera and so align the photometry using 
a large number of relatively bright stars. If the photometric error were reduced from the 
present 5.5% to 2%, then the mass error would be reduced from 17% to 15%. If this were 
combined with a 50 /xas parallax measurement, the mass error would be reduced from the 
above-mentioned 10% to 7%, while perfect knowledge of the distance would, under these 
circumstances, reduce the mass error to 4%. 

The third path would be improved microlensing data reductions. All the microlensing 
analysis has been conducted on the basis of the original MACHO SoDoPHOT pipeline pho- 
tometry. This pipeline produced 10 10 measurements of very high quality, but with modern 
image-subtraction routines, it may be possible to do better. However, since the SoDoPHOT 
errors have been renormalized (see § 3.1) some of this improvement has already been achieved. 

A.G. was supported by grant AST 02-01266 from the NSF. D.P.B was supported by 
grants AST-0206189 from the NSF and NAG5- 13042 from NASA. 
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Fig. 1. — Likelihood contours in the 7r E plane shown at A% 2 = 1, 4, 9, 16, 25, 36, and 49 
relative to the minimum. This should be compared to Figure 3 of Gould (2004). The errors 
here are smaller, partly because of error renormalization and partly because of the addition 
of HST baseline photometry of the source from Alcock et al. (2001). 
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Fig. 2. — Likelihood contours (A% 2 = 1, 4, 9) for the log mass M and distance modulus (M — 
m)o of MACHO-LMC-5 based on the lightcurve of the event, the source- flux measurement 
of Alcock et al. (2001), and constrained by the proper-motion measurement of Drake et al. 
(2004). Each set of contours is shown relative to its own minimum. The left-hand minimum 
is actually higher than the one at the right by A% 2 = 6. The vertical lines show the best-fit 
distance modulus and 1 o confidence interval derived from the trig parallax measurements of 
Drake et al. (2004), while the point with error bars shows logM and (M — m)o as determined 
photometrically from HST WFPC2 data by Alcock et al. (2001). Note that the photometric 
(contours) and astrometric (vertical lines) determinations of (M — m) are in agreement at 
better than 1 a. 



-20- 




Fig. 3. — Similar to Fig. 2 except that the microlensing/proper- motion determination has 
now been combined with the lens-source relative parallax measurement of Drake et al. (2004). 
The errors are 17% in the mass and 6% in the distance. There is good agreement with the 
mass and distance estimates based on HST photometry, shown as a point with error bars. 



